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Abstract. In this paper we construct the main algebraic and differential properties 
and the weight functions of orthogonal polynomial solutions of bivariate second-order 
linear partial differential equations, which are admissible potentially self-adjoint and of 
hypergeometric type. General formulae for all these properties are obtained explicitly 
in terms of the polynomial coefficients of the partial differential equation, using vector 
matrix notation. Moreover, Rodrigues representations for the polynomial eigensolu- 
tions and for their partial derivatives of any order are given. Finally, as illustration, 
these results are applied to specific Appell and Koornwinder orthogonal polynomials, 
solutions of the same partial differential equation. 
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I. Introduction 



The theory of orthogonal polynomials in one variable is in permanent expansion due 
to its relationship with other areas of mathematics and also with several applications 
in physics and engineering. They provide a natural way to solve many types of impor- 
tant differential equations of mathematical physics, expanding solutions in appropriate 
Fourier series of orthogonal polynomial basis. They play therefore an important role 
^ ■ in the study of wave mechanics, heat conduction, electromagnetic theory, quantum me- 
chanics or mathematical statistics. 

In this context, it is first relevant to study whether a (one variable) polynomial family, 
{p n (x)} neN (x G E), is orthogonal. This problem is solved in different ways but the 
Favard's theorem [20], linking orthogonality and the fundamental three-term recurrence 
relation 



(1) xp n (x) = a n p n+1 (x) + (5 n p n (x) + 7„p n _i(x) , 7„ ^ 0, 

provides certainly the most powerful characterization [7]. Here, if p n (x) = g n>n x n + 

9n,n-ix n ~ l + g n ,n-2X n ~ 2 H , then orthogonality of the p n -family leads easily to the well 

known expressions (see e.g. [7j, [U], [2U equation (1.4.17), p. 14] or [251 P- 36]): 

/r,\ 9n,n 9n,n—l &n9n+l,n 9n,n—2 Pn9n,n—1 (%n9n+l,n—l 

(2) a n = , fi n = — '- , 7„ = — : '■ '■ . 

9n+l,n+l 9n,n flVi— l,n— 1 

It is also important to provide ways of constructing efficiently these polynomials. For, 
several approaches are at hand. Besides the use of the recurrence (CQ) itself, we can 
mention (in a non-exhaustive way) the generating function methods, those based on a 
Rodrigues formula or the hypergeometric approach which gives nice and useful represen- 
tations of the polynomials in terms of hypergeometric series. 
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As it is very well known, among all the one variable orthogonal polynomials, the four 
classical continuous families of Jacobi, Laguerre, Hermite and Bessel, are those sharing 
the widest set of properties. Besides the three-term recurrence (JTJ) [TIESIED], they can be 
characterized in a number of ways, e.g. they are orthogonal polynomial solutions of the 
hypergeometric type differential equation [6j [25] and the fc-th derivatives of each family 
are again orthogonal and belong to the same family [TJ [25]. Moreover, the orthogonality 
weight functions satisfy Pearson-type equations [6], [21] giving rise to Rodrigues formulae 
[H [25] for the corresponding orthogonal polynomials and for their derivatives of any 
order. Also, the orthogonal polynomials satisfy a number of algebraic and differential 
properties such as derivative representations [H [23] or also structure relations [TJ [T7] , 
among other properties. The list of references in this paragraph is not exhaustive but 
only indicative of the kind of references that could be examined on this topic. 

In these classical settings, it is remarkable that the coefficients appearing in all the 
aforementioned algebraic and differential characterizations can be explicitly computed in 
terms of the polynomial coefficients cr(x) and t(x) of the hypergeometric-type differential 



Our main contribution is to extend this remarkable property to the bivariate situation 
for polynomial solutions of admissible potentially self-adjoint partial differential equa- 
tions of hypergeometric-type [201 I2H [22]. These polynomials play an important role in 
many applications, for instance in spectral//ip-finite element methods for solving partial 
differential equations [H [13] . 

One essential difference between polynomials in one variable and in several variables 
is the lack of an obvious basis in the latter [9]. One possibility to avoid this problem is 
to consider graded lexicographical order and use the matrix vector representation, first 
introduced by Kowalski [181 HE] an d afterwards considered by Xu [311 [35]. I n i&ct, using 
this point of view, in [2] the authors proved some structure and orthogonality relations for 
the successive partial derivatives of the vector orthogonal polynomials associated with 
a quasi-definite moment functional which satisfies a Pearson-type partial differential 
equation. 

In this paper we deal with bivariate polynomials written in vector representation (and 
graded lexicographical order) which are solutions of admissible potentially self-adjoint 
linear second order partial differential equation of hypergeometric type. In this context, 
we prove that (as it happens in the one variable hypergeometric type case) the coefficients 
characterizing the three-term recurrence relations, the first structure relations and the 
derivative representations fulfilled by the vector polynomials can be written explicitly in 
terms of the coefficients of the partial differential equation they satisfy. In the bivariate 
discrete case some results in this direction have been already given in [261 [23 [2H] ■ 

The structure of the paper is as follows: In Section [21 after introducing basic defini- 
tions and notations, we present in Proposition I2.5l the general framework to be considered 
through the paper, i.e. admissible potentially self-adjoint second order partial differential 
equation of hypergeometric type. In Section [3] we give the partial differential equations 
for the partial derivatives of the eigensolutions and we construct the corresponding weight 
functions for the orthogonal polynomials. Then, the relations linking these weight func- 
tions are obtained and they allow us to deduce a Rodrigues formula for the orthogonal 
polynomial solutions and for their partial derivatives of any order. In Sections H] and 



equation Q21 [23 (23 [361 EH [38. 
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using vector matrix notation [91 US [19] , general formulae for the main algebraic and 
differential properties (three-term recurrence relations, structure relations and derivative 
representations) are explicitly obtained in terms of the coefficients fully characterizing 
the partial differential equation, both in monic and non-monic cases. Finally, in Section 
[6] our results are applied to three different eigensolutions of the same partial differential 
equation, which are orthogonal on the same domain with respect to the same weight 
function g(x) = x^ 1 y 13 ^ 1 and so, they form biorthogonal families. First, the two pa- 
rameter monic Appell polynomials family [11] is analyzed in detail. From that, for the 
two remaining cases (non-monic Appell polynomials [I] and Koornwinder polynomials 
[T5] ) the corresponding properties can be deduced in two ways: from the results in Sec- 
tion H] or by means of their expansions in terms of monic Appell polynomials which we 
establish here. 

2. Vector representation and admissible partial differential equations 

OF HYPERGEOMETRIC type 

Let x = (x, y) G M 2 , and let x n (n G No) denote the column vector of the monomials 
x n ~ k y k , whose elements are arranged in graded lexicographical order (see [91 p. 32]): 

(3) x n = (x n ~ k y k ) , < k < n, n G N . 

Let {Pn-k fc( x ' v)} be a sequence of polynomials in the space IT 2 of all polynomials of total 
degree at most n in two variables, x = (x,y), with real coefficients. Such polynomials 
are finite sums of terms of the form ax n ~ k y k , where a G R. 

From now on the (column) vector representation |T8l US] will be adopted, so that P n 
will denote the (column) polynomial vector 

(4) P n = (P n y*, y), PZ_ hl (x, y),..., P£ n _i(*, y), P o yz, y)) T 
Then, each polynomial vector P n can be written in terms of the basis ([3]) as: 

(5) P n = G ntn X n + G n<n -iX n 1 + • ■ ■ + G n fl X°, 

where G n j are matrices of size (n + 1) x (j + 1) and the leading matrix coefficient G n ^ n 
is a nonsingular square matrix of size (n + 1) x (n + 1). 

Definition 2.1 (Monic polynomial vector). A polynomial vector P n is said to be monic 
if its leading matrix coefficient G n>n is the identity matrix (of size (n + 1) x (n + 1)); i.e.: 

(6) P n = x" + G^.iX"" 1 + • • ■ + G nfi x° . 
Then, each of its polynomial entries P™_ k k (x,y) are of the form: 

(7) Pn-k fc( x ) y) = x n ~ k y k + terms of lower total degree . 
In what follows the "hat" notation P n will represent monic polynomials. 

Definition 2.2 (Orthogonality). Let £ be a moment linear functional acting on II 2 . A 
sequence of polynomials {P^-k k( x ^v)} C LT 2 {n G N ), is said to be orthogonal with 
respect to C or, equivalently, {P n }n>o (as defined by Eqs. (|4])-([5])) is a vector orthogonal 
polynomial family with respect to C, if for each n G No there exist an invertible matrix 
H n of size n + 1 such that: 

(8) £ [(x m PD] = G M {m+1 ' n+1) , n > m, 

(9) C [(x n P^)] =H n e M {n+1 ' n+1) . 
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An integral representation of this orthogonality functional £ can be constructed by 
means of a weight function g := g(x, y) defined in a certain domain Del 2 : 



which is defined in the set II 2 provided that all the above integrals exist. Then, the 
family {P n } n >o is said to be orthogonal with respecto to g in the domain D. 

In this multivariate context, Bochner [6] posed the problem of identifying those fami- 
lies of polynomials which are eigenf unctions of a second order linear partial differential 
operator. Krall and Sheffer [2D] started to study eigenf unctions which are orthogonal 
over a domain giving conditions of admissibility and a first attempt of classifying admis- 
sible equations. Engelis [10] gave a detailed list of second order linear partial differential 
equations for which orthogonal polynomial in two variables are solutions. This question 
was afterwards studied and systematically described by Suetin [32] . In this paper, we an- 
alyze admissible potentially self-adjoint partial differential equations of hypergeometric 
type. 

In order to present this study, we consider a bivariate class of linear partial differential 
equations, introduced as "the basic class" by Lyskova [22] for the multivariate case (see 
also [3]) and called here hypergeometric type equations: 

(11) (a x x 2 + hx + c x )d xx u(x, y) + 2(a 3 xy + b 3 x + c 3 y + d 3 )d xy u(x, y) 

+ (a 2 y 2 + b 2 y + c 2 )d yy u(x, y) + (etx + fi)d x u(x, y) + (e 2 y + f 2 )d y u(x, y) + Xu(x, y) = 0, 

where aj, bj, Cj, dj and A are real numbers. The solutions of this equation have the 
remarkable property that all the partial derivatives of any order of these solutions are 
also solutions of an equation of the same form. 

Moreover, we shall also consider admissible partial differential equations. 

Definition 2.3. A second order partial differential equation is admissible if and only if 
|20[ 132] for any non-negative integer n there exists a number A„ such that equation ffTTT) 
with A := A„ has n + 1 linearly independent solutions which are polynomials of total 
degree n and has no non-trivial solutions in the set of polynomials of total degree less 
than n. 

The following characterization has been proved in [221 E2] • 

Proposition 2.4. Equation 177]) is an admissible second order partial differential equa- 
tion of hypergeometric type if and only if it can be written in the form 



(ax 2 + hx + ci)d xx u(x, y) + 2(axy + b 3 x + c 3 y + d 3 )d xy u(x, y) 

+ (ay 2 + b 2 y + c 2 )d yy u(x, y) + (ex + f!)d x u(x, y) + (ey + f 2 )d y u(x, y) + X n u(x, y) = 0, 

where X n = —n((n — l)a + e) and the coefficients a, bj, Cj, d 3 , e, fj are arbitrary fixed real 
numbers, but the numbers a and e are such that the condition 



(10) 




(12) 




Wk '■= ak + e 7^ 0. 



In these conditions an orthogonality weight function defined in certain domain of 
IR 2 which is related with the partial differential equation (fT2"|) can be given for their 
orthogonal polynomial solutions as shown in [32J. 
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Proposition 2.5. Let a(x,y) be the discriminant of Eq. ( LU| ) ; i.e.: 

(14) a(x, y) = (ci + x (b\ + ax)) (c 2 + y (b 2 + a y)) - (d 3 + b 3 x + (c 3 + a x) y) 2 , 
and DcK 2 be the domain: 

(15) D = {(x, W )eM 2 : a(x } y)^0} . 
Defining the two functions [321 eq. (15), p. 132] 

(16) /3(x, y) = (-&i - c 3 + ft - 3ax + ex)(ay 2 + b 2 y + c 2 ) 

- {-b 2 - b 3 + f 2 - Say + ey)(axy + b 3 x + c 3 y + d 3 ) , 

(17) 7(x, y) = -(ci + x(6i + ax))(& 2 + &3 - /2 + 3ay - ey) 

+ (6i + c 3 - /i + 3aa; - ex)(d 3 + 6 3 x + (c 3 + ax)y) , 
assuming that in D the following condition holds true: 
/io\ 5 / 7(^,Z/)\ <9 / /3(a;,2/) > 



(9a; \a(x,y) J dy \a(x,y) 

considering the weight function or integrating factor of (W) given by [32], Equation (22), 
P- 134] 



(19) g(x,y) = exp{ [ 1 ^' y dy + 



P{x,y) 

a(x,y) 



y=yo 



dx 



which determines (up to a multiplicative constant) the functional 
(20) C(P) = JJ P(x,y) e (x,y)dxdy, Pen 2 , 

defined in the set provided that all such integrals exist, then, there exists a unique 
monic vector polynomial family {P n }n>o solution of |Hj) and orthogonal with respect to 
q in D, i.e. satisfying 



e .M (m+1 ' n+1) , ifn>m, 



(21) // x m FT Q (x,y)dxdy 
l> 1 H n e M {n+1 ' n+l) , if m = n, 

where H n {of size (n + 1) x (n + 1)) is nonsingular. 

Remark 1. From the partial differential equation (fl2"|) it is possible to introduce the linear 
operator D by 

(22) Vu = (ax 2 + b\x + c\)d xx u + 1(axy + b 3 x + c 3 y + d 3 )d xy u 

+ (ay 2 + b 2 y + c 2 )d yy u + (ex + fi)d x u + (ey + f 2 )d y u, 

which is potentially self-adjoint if and only if (IT51) holds true [32, Theorem 1, p. 133]. 

So, as it has been mentioned in the introduction, in this paper we shall be concerned 
with bivariate orthogonal polynomial families {P n } n >o (in the sense of Definition 12.21) 
written in vector representation dH)-©, which are solutions of admisible potentially self- 
adjoint second order partial differential equations of hypergeometric type characterized 
in the Proposition 12.41 
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3. Weight functions and Rodrigues formula for the polynomials and 

for their partial derivatives 

By differentiating (TT2"|) (r + s) times, it turns out that 

^ S) (x,y) = ^(x,y), r,s = 0,1,2,..., 
satisfies an admissible second order partial differential equation of hypergeometric type 

(23) (ax 2 + b& + c 1 )d xx z ( - r ' s \x,y) + 2(axy + b 3 x + c 3 y + d 3 )d xy z {ryS) (x, y) 

+ {ay 2 + b 2 y + c 2 )d yy z^ s \x, y) + r x r ' s \x)d x z^ s \x, y) 

+ (y)d y z^ (x, y) + fi r+a z™ (x, y)=0, 

where 

(24) 4 r ' s) (x) = (e + 2a(r + s))x + fi + rb x + 2sc 3 , 

(25) r^ s \y) = (e + 2a(r + s))y + f 2 + 2rb 3 + sb 2 , 

(26) [i r+s = X n + (r + s)e + (r + s)(r + s - l)a, 

and d 3 , e, /j are arbitrary fixed real numbers satisfying conditions ([TBI . Observe 

that equation (1231 and the above relations have been obtained in the multivariate case 
for hypergeometric type equations (not necessarily admissible) in [221 Theorem 1]. 
Let us also introduce the linear operator 

(27) V^'^z = (ax 2 + b x x + ci)d xx z + 2(axy + b 3 x + c 3 y + d 3 )d xy z 

+ (ay 2 + b 2 y + c 2 )d yy z + r^ 3 \x)d x z + r^ s \y)d y z. 

Applying condition flJSJ) to the operator (j2"7j) . we obtain that this operator is potentially 
self-adjoint in a domain D, if and only if 



(28) 


d ( u(x,y)\ 
dx \a(x, y) ) 


\+(s- 


- r) TTx{ 


1 da(x,y)\ d 


a(x,y) dy J dy 


where 










(29) 






a(x,y) 


= AC-B 2 , 


(30) 






u(x,y) 


ox ox 


(31) 






9(x,y) 


= 2C— -5— , 

dy dy ' 


and 










(32) 




A = 


A(x,y) 


= ax 2 + b\x + ci, 


(33) 




B = 


B(x,y) 


= axy + b 3 x + c 3 y + d 3 , 


(34) 




C = 


C(x,y) 


= ay 2 + b 2 y + c 2 . 



Q(x,y) 

dy \a(x,y) 



0. 



If the operator T)^'^ defined in ( J27j) is potentially self-adjoint in a domain D, there 
exists in this domain a positive and twice continuously different iable function g^'^ (x, y) 



BIVARIATE PARTIAL DIFFERENTIAL EQUATIONS AND OP SOLUTIONS 

which is the solution of the system of differential equations (Pearson type equations) 
equations (7) and (8)] and [321 p. 132] 

1 dg( r < s \x,y) ft r > s \x,y) 



(35) 



g( r,s ^(x,y) dx a(x,y) 
1 dg( r > s \x,y) ^ r ' s \x,y) 

g( r > s )(x,y) dy a(x,y) 



where 

(36) P {r > s \x,y) = I3(x,y) + r-£(x,y) + s6(x,y), 

dot 

(37) 7 (r ' s) (x, y) = 7(2, y) + ru(x, y) + s— (x, y), 

dy 

and the polynomials (3(x,y) and j(x, y) have been defined in ( fl6l) and (|T7j) respectively 
From the Pearson type equations (133]) we obtain 



(38) g ir ' s \x,y) = exp{ f ^—^f}- dy + 



y ^-\^jyj Jx 



(3 {r ' s) (x,y) 
a(x,y) 



y=yo 



dx 



up to a multiplicative constant. 

Now, we are in a position to establish the relation linking g( r,s \x,y) and g(x,y). 

3.1. Relation between weight functions. Now, we can establish the connection be- 
tween y) and g^ fi \x,y) = g(x,y), given in (138]) and ( 119]) respectively. From 
equation ( 138|) . after straightforward computations, we obtain that 

(39) g^ s \x, y) = 0^(x, y)g(x, y), r, s = 0, 1, 2, . . . , 

up to a multiplicative constant, where (j)^ r ' s '(x, y) is a polynomial which explicit expression 
depends on the coefficients of the partial differential equation ([12]) . After solving the non 
linear system of equations (128]) for any r and s, we can reduce the solutions of the system 
to the following ten cases: 

(i) If bi = 2c3 and b 2 = 2b 3 , we have 

<f>( r > s \x,y) = [a(x,y)] r+s , 

where 

a(x,y) = -(d 3 + b 3 x + (c 3 + ax)y) 2 + (c x + x(2c 3 + ax))(c 2 + y(2b 3 + ay)). 

(ii) If c 3 ^ 0, d 3 ± 0, 6 3 + 0, a = , d = ^f^, and c 2 = ^J^, we have 

* ( //j (d 3 + C 3 ^(d3 + M S ' 

where 

03C3CX3 

x ((4 + 63^X^3 + c 3 y){-b x {^ 2 d 3 - b 3 d 3 + 6 3 c 3 y) + c 3 {b 2 {d 3 - b 3 x) + 26 3 (6 3 x + c 3 y)))). 

(iii) If a = bi = c\ = c 3 = 0, we obtain 

(r ' s) (x, y) = [a(x, y)] r , a{x, y) = (d 3 + b 3 x) 2 . 
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(iv) If a = b 2 = b 3 = c 2 = 0, we have 

(f)^ s \x,y) = [a(x,y)} s , a(x,y) = (d 3 + c 3 y) 2 . 

(v) If a = b 3 = c 3 = d 3 = 0, we have 

^ y) = (ci+ '^7w a(x ' v) = (C1 + blX)[a - + biyY 

(vi) If a 7^ 0, 63 = c 2 = d 3 = 0, and C\ = ( - bl ~^° 3 ; we obtain 

V ; y r (c 3 + ax)*' 

where 

_ (c 3 + ax)y(b 2 {bi - c 3 + ax) + a(6i - 2c 3 )y) 

a 

(vii) If C3 7^ 0, a = 63 = 0, 61 = C3, and c 2 = we obtain 
where 

/ \ (4 + c 3 y)(& 2 (ci + c 3 x) - c 3 {d 3 + c 3 y)) 

a(x,y) = . 

c 3 

(viii) If 63 7^ 0, a = c 3 = 0, 6 2 = b 3 , and c\ = bid 3 /b 3l we obtain 

]r+s 



)= M*,v)Y 

(d 3 + b 3 xY 



where 

(d 3 + b 3 x)(-b 3 (d 3 + 6 3 )x + 6 x (c 2 + 63?/)) 



a(x,i/) 



(ix) If a 7^ 0, ci = c 3 = d 3 = 0, and c 2 = 12 , we obtain 

]7~ + S 



where 

_ x(b 3 + ay){a{b 2 - 2b 3 )x + 61 (6 2 - 63 + ay)) 

a 

(x) If ci = c 2 = d 3 = 63 = c 3 = 0, we obtain 

lr+s 



x s y r 
where 

y) = xy{ab 2 x + 61(62 + ay)). 

Remark 2. Observe that in all the above cases the polynomial (f)( r,s \x,y) defined in fl39l) 
can be factorized as 

(40) ^ s \x,y) = [<f>^(x,y)Y[^ 1 \x,y)}\ 

for any r and s. 
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Remark 3. It is important to notice here that different cases could give rise to the 
same (p^ r ' s \x,y) associated with the same partial differential equation. This situation 
appears in the example studied in detail in Section [6] where the coefficients of the partial 
differential equation ( 1541 fulfill the conditions of cases (vi), (ix) and (x). 

Remark 4. We must now mention that in [22] Theorem 3], Lyskova presented the equa- 
tion f[3"9"]) in a non-explicit form. 

3.2. Rodrigues formula. One of the main problems in the theory of orthogonal polyno- 
mials in several variables is to obtain explicit expressions for the orthogonal polynomial 
solutions of the partial differential equation. In this direction we could mention the 
works of Engelis [10] , who derived the Rodrigues formula for some classes of orthogonal 
polynomials in two variables, and Suetin [32] Theorem 3, p. 151], who showed that 
this Rodrigues representation is one of the ways of constructing explicitly orthogonal 
polynomial families in the the potentially self-adjoint case. 

As a consequence, in the context considered here (admissible potentially self-adjoint 
partial differential equations of hypergeometric type), we have the following explicit 
Rodrigues formula for the polynomial solutions of (I12p of total degree n + m 

(x,y)[^\x,y)} n [4 > ^(x,y)} m , 

where K„ im is a normalizing constant and the polynomials (j)( r ' s \x,y) have been intro- 
duced in (JSnD and (jiO]). 

Having in mind that the partial differential equation for the partial derivatives (1231 
is also of hypergeometric type with the same discriminant a(x,y) defined in (1T4"]) . we 
also have a Rodrigues representation for the partial derivatives of any order, polynomial 
solutions of (1231 . given by 
(42) 

where g^ r ' s \x,y) is given in (138]) and N n ,m,r,s is a normalizing constant. In this way we 
have obtained a natural extension to the bivariate case of the Rodrigues representation 
for classical orthogonal polynomials in one variable [24"] . 



(41) 



O ' 



nix. 11) dx n dii r 



n,m,r,s 



o( r >*)(x.v) dx n - r dv r 



4. Explicit expressions for algebraic properties 

Let us consider a vector polynomial family {P n } ng N solution of (1T21 orthogonal in 
the sense of Proposition 12.51 i.e. it is orthogonal with respect to a weight ( 1191) and 
satisfies (12ip in an appropriate domain D G R 2 . In such a conditions it can be proved 
that the family satisfy a number of algebraic and differential properties. Here we focus 
our attention in three of the most relevant: the three-term recurrence relations, the 
structure relations and the derivative representations. As already mentioned, our aim 
in this section is to give all of these relations in terms of the matrix coefficients in the 
expansions (J5J) of the {P n }- family elements. 

For, let us first introduce the matrices L n j of size (n + 1) x (n + 2) defined by 

(43) Ln,ix n+1 =«", L n , 2 x m+1 = y x n , 
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(I 



L n ,i 



o °\ 



and L 



n,2 



/0 1 



V O 1 °/ 

Let us observe that 
(45) x 2 x" = L nil L n+lil x"+ 2 , y 2 x" = L ni2 L n+1 , 2 x"+ 2 



1/ 



and for j = 1, 2, 
(46) 



Ln,2Ln+l,l 



T ■ T/ 



L n ,\L 



n+1,2, 



+1) 



where denotes the identity matrix of size n + 1. 
Moreover, for n > 1, 



(47) 



cU n = E nil x"-\ 

a / x n = E n , 2 x n - 1 , 



where the matrices E n j of size (n + 1) x n are given by 

/On o ^ 

n- 1 



(4J 



O 
\0 ... 



and E n 2 



0/ 



/o 
1 



V O 



o\ 



O 



4.1. The three— term recurrence relations. The existence of a recurrence relation 
for a vector orthogonal polynomial family can be established in more general settings 
than those considered here. More precisely, the following existence theorem is proved in 

m- 

Theorem 4.1. Let C be the positive definite moment functional as defined in ^^j and 
{En}n>o be an orthogonal family with respect to C Then, for n > 0, there exist unique 
matrices A n j of size (n + 1) x (n + 2), B n j of size (n + 1) x (n + I), and C n j of size 
{n + 1) x n, such that 



(49) 



^nj'IPn+l + B n jW n + C nJ 'P n _i, j — 1,2 



with the initial conditions P_i = and Pq = 1. i/ere i/ie notation x± = x, x 2 = y is 
used. 

Now it is possible to generalize the well-known expressions (T5]) for the one variable 
case to the bivariate case. This is done in the following proposition which is proved with 
the help of the auxiliary matrices L n j defined in (j43p — (j44p . 

Theorem 4.2. The explicit expressions of the matrices A n j, B n j and C n j (j = 1,2) 
appearing in UM in terms of the values of the leading coefficients G n>n , G njn _i and 
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G n ,n-2 in the expansions |3P are given by 

A n j = Gn,nL n jG n+1 ^ n+ i, n > 0, 
Bqj = —A jG\ : o, 

(50) < B n j = (G n ^iL n _i j — A n jG n+ i tn )G n ^ n > 1, 

G\j = — (Ax j(j2,o + BijGifi), 
^G n j = (G n>n -2L n -2,j — A n jG n+ i tn -i — B n) jG n)n —i)G n _i >n _i, n > 2 . 

Proof. In equation ( )49l) . substitute P n as given in ((5]), equate the coefficients of x fc for 
k = n, n — 1, n — 2 and solve the corresponding linear system. □ 

The above result is valid for any orthogonal polynomial sequence (12 ip . From now on 
we shall consider algebraic properties involving the partial derivatives of the orthogonal 
polynomials. To set up sufficient conditions for the existence of this type of relations is 
not enough with the orthogonality of the vector polynomial family. For these relations 
to exist, the orthogonal polynomial family has to be solution of an admissible potentially 
self-adjoint second order partial differential equation of hypergeometric type (JI2]). With 
these assumptions, next we present explicit expressions for the three-term recurrence 
relation of the first partial derivatives, structure relations and derivative representations. 

4.2. The three— term recurrence relations for the first partial derivatives. Let 

us define, for j — 1, 2 

(51) QW = L n ,^ n+1 = (^Kofo tf), ^.y), • • • , £.p£\x,y)) T 

where P n is defined in (@J-([5]) and L n j are defined in (jUJ. Here the notation x\ = x, 
x 2 = y is used. Thus, we have 

(52) ^ Pn+1 = Grf S x ™ + G ^n-i* n ~ l + ■■■ + Gd% x°, 
where 

(53) Gd^ k = G n+1<k+1 E k+1J , 0<k<n, j = 1,2, 

are matrices of size (n + 2) x (k + 1). 

Let L be the positive definite moment functional as defined in ( 120]) and {P n }n>o be an 
orthogonal family with respect to C, so that it satisfies the admissible partial differential 
equation of hypergeometric type (fl2l) . 

Then, for n > and j = 1,2, there exist matrices A^j of size (n + 2) x (n + 3), B^'- 
of size (n + 2) x (n + 2), and C^j of size (n + 2) x (n + 1), such that 

(54) .o-^-i -^ ? ^ • B S^- p ^ • r ^jj7 j ? - i = 1 > 2 > 

with the initial conditions ^§-Po — and §|-Pi = Gi,iEij. 

The explicit expressions of the matrices A^, B^l and C^j can be obtained in terms 

of the values of the leading coefficients Gdn} n , Gd^ n _x and Gd^' n _ 2 , defined in (|53|) . using 
the Proposition 12.51 but they are not unique. 
Since 

(55) ^-P n+ i = Ll0\ n>0, 
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for j = 1, 2 we can write a three-term recurrence relation for the Qn polynomials, which 
is now unique after multiplying by the • matrices. In this way, the following theorem 
can be proved. 



Theorem 4.3. Let C be the positive definite moment functional as defined in §20y and 
{Pn}n>o be an orthogonal family with respect to C Then, for n > and j = 1,2, there 
exist unique matrices A^- of size (n + 1) x (n + 2), B^j of size (n + 1) x (n + I), and 
C^l of size (n + 1) x n, such that 



(56) 




-B U) Q {j) +C (J) Q U \, j = 1,2 


where 






(57) 
(58) 
(59) 


50') 

B n,j 


- T -A {j) T T 

T r0) r r 

- - Lj n,j J -> n j- L ^n t j, 

- r rO) T t 



with the initial conditions Q_\ = (0) and = L jGdi^. Here the notation x\ = x, 
X2 = y is used. 

4.3. First structure relations. In the one variable case, the so-called first structure 
relation plays an important role since, e.g. it gives rise to lowering and rising operators. 
Next we present the extension to the bivariate case. 

Theorem 4.4. Let {P n }neN be a vector orthogonal polynomial family satisfying Propo- 
sition 001 Then, for n > 1, there exist unique matrices W n j of size (n + 1) x (n + 2), 
S n ,j of size (n + 1) X (n + l),and T nj - of size {n + 1) x n, such that 

(60) 0,-O, V)^n = W n:J F n+1 + S nJ F n + T nij P n _!, j = 1, 2, 
where 

(61) fa(x,y) = (1 ' o) (:r,y) , 2 (x,y) = ^ , 

and i/ie polynomials 4>^ r,s \x,y) have been introduced in and (f^^- Wit/i i/je notation 

(62) <l>j(x,y) = ajX 2 + fyxy + <yjy 2 + 5jX + Ejy + Uj, j = l,2, 
i/ie coefficients of the structure relations are explicitly given by 

,n+l Q ' n,n^n. j {ctjL n -i,iL n i + j3jL n _i^L n ^ + 7j-^n-i, 2-^71,2) j 

+ C n) „_iE„_ij (ajL n _2,i-L n _i ) i + /3jL„_2,2-^n-i,i + 2-^-1,2) , 

TnjG n — X,n— 1 UjGn^nEnj W n jG n -\-i^ n — \ S n jG ntn —i 
+ Gn^-iE^.ij {b~jL n _2 : \ + £jL n _2,2) 

2^n-2,j (c^jL n _ 31 L n _2 t l + /3jL n _ 3:2 L n _2 : l + Jj -^n-3,2-^n-2,2) • 

Proof. Since 
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are polynomials of total degree n + 1 we can write for j = 1,2, 

d n+1 
(f> j {x,y)—¥ n = Y^KjPt, n>l. 

j e=o 

As a consequence of Section El the polynomials -§^n are orthogonal with respect to 
<fii(x,y)g(x,y), and are orthogonal with respect to (p 2 (x,y)g(x,y). If we multiply 
from the right the latter equation by P^ and apply the corresponding orthogonality we 
get A Uj £ = for £ < n — 1. 

The explicit expressions for the matrix coefficients are obtained from equation (|60|) . 
by substituting P n as given in (jSJ), equating the coefficients of x fc for k = n, n — 1, n — 2 
and solving the corresponding linear system. □ 

Observe that if we know the explicit expression of the polynomials <pj(x,y), j = 1,2, 
by using the results given in Section 13. 1[ it is possible to obtain explicitly the matrices 
W n j, S n j and T n j in fl60|) in the following way: substitute P n as given in (jSJ) equate the 
coefficients of x fc for k — n, n — 1, n — 2 and solve the corresponding linear system, by 
using (jl3D, (USD and (HH). 



4.4. Derivative representations (or second structure relations). Next we present 
finite-type relations between the orthogonal polynomial sequence {P n } and the sequence 
of the partial derivatives {^-P n }. 

Theorem 4.5. Let {P n }neN be a vector orthogonal polynomial family satisfying Propo- 
sition \2.5[ Then, for n > 2 we have 

d d d 

(63) P n = V n j— — Pn+i + Y n j— — P n + Z n j— — P n -i, 

dxj dxj dxj 

where the matrices V n j of size (n+1) x (n + 2), Y n j of size (n + 1) x (n + 1), and Z n j 
of size (n + 1) x n are given by 

V ■ — A ■ - 

Y — B ■ - 

and the matrices A n j, B n j and C n j are given in £5U\) and the matrices A^_ X a, B^_ x - 
and C„_i j are introduced in |5^| ). 

Proof. The above result is a consequence of ( )49|) and (|54|) . □ 
Moreover, from (155|) we can deduce the following relation: 

(64) P^K.^QW + ^i^^J+^i^.QW, n>l, j = l,2, 
where L£ • is the transpose matrix of L n j defined in 
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5. Explicit expressions for algebraic properties in the monic case 



As pointed out in the introduction, in this section we give the explicit expression 
of the coefficients G njn _i and G„ in _ 2 in fl6]) in terms of the coefficients a, bj, Cj, d 3 , e, 
fj fully characterizing the already mentioned partial differential equation f fT2|) . After 
that, results already given in the previous section will allow us to express (for monic 
polynomials) the three algebraic and differential properties here considered ( 1491) . ( jBTJl) 
and f lBHj) . in terms of the admissible partial differential equation coefficients in ffT2"j) . 



Proposition 5.1. Let P n (n G No) be a monic vector polynomial, as given by the ex- 
pansion (TJJ), solution of an admissible hypergeometric type partial differential equation of 
the form (TUP- Then, the matrix coefficients G n , n _i G Al (n+1 ' n) and G n , n _ 2 G A^ (n+1 ' n " 1) 
in (0|) can be written in term of the coefficients d 3 , e, fj in ([X|) as: 



92,1 92,2 



(65) 



v o 



9n-\,n-2 9n-\,n-\ 

9n,n—l 9n,n 
9n+l,n/ 



(n > 1) 



where, for 1 < i < n, 



9i,i 

9i+l,i 



( n + i_,)(( n _,) &1 + 2 (^-i)c 3 + / 1 ; 

^2n-2 

i((i-l)b 2 + 2(n-i)b 3 + f 2 ) 

^2n-2 



and 



(66) 



/ 9i,i 

92,1 
93,1 



G 



n,n—2 



o \ 



92,2 
93,2 



93,3 



9n-l,n-3 9n-l,n-2 9n-l,n-l 
9n,n—2 9n,n—l 
9n+l,n-lJ 



(n > 2) 
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where, for 1 < i < n — 1, 

(n — i)(n + 1 — i) 



9i 



x (n7 2n _ 2 d + ({n - i)h + 2(i - l)c 3 + /i)((n - z - l)&i + 2(i - l)c 3 + A)), 



Z ( Tl — Z ) 

= K - J —(fif 2 + + h{2{n - 2 + 2{i - 2){n - i - l))c 3 

^2n-2^2n-3 

+ (2n - 2% - \)h) + (2i - l)c 3 / 2 + (i - l)fo 2 ((2z - l)c 3 + /i) 
+ („ _ i _ _ i) &2 + ( 2 „ - 2i - 1)6 3 + / 2 )), 

z(i+l) 
5 , j+2,i — ~ 

^^2n-2^2n-3 

x (w 2 n- 2 C2 + ((i - 1)6 2 + 2(n - i - 1)& 3 + / 2 )(z&2 + 2(n - i - 1)6 3 + / 2 )), 
In a// 0/ t/iese expressions vj n = na + e 7^ as already shown in Eq. [W\l . 

Proof. Plug into equation f[T2"j) the expansion (jH]) and then make equal zero the coefficients 
of the column vector of monomials x fc (defined in (jHJ)) for = n, n — 1, 71 — 2. □ 

Now, having in mind that in the monic case G n ^ n = I n +i, from these results and 
Theorem I4.2[ we can deduce the following corollaries. 

Corollary 5.2 (Three-term recurrence relations). For monic polynomials and j = 1,2, 
the coefficients of the three-term recurrence relation fijtfy are given in terms of the coef- 
ficients of the second order partial differential equation by 

(67) A n j = L n j, 

(68) 5b,i = , £0,2 = {-^ 

(69) B n j = G n}n -iL n -ij — L n jG n+ i^ n , n > 1, 

t -cie 2 +/i(bie-a/i) \ / -d-j ff+b 3 ef 1 +c 3 ef 2 -af 1 f 2 * 

e 2 (a+e) \ _ I e 2 (a+e) 

-d 3 e 2 +b 3 ef 1 +c i ef2-af 1 f2 J ' 0l - 2 I -c 2 e 2 +f 2 (b 2 e-af 2 ) 
e 2 (a+e) J \ e 2 {a+e) 

(71) C n j = G J n>n _ 2 X n _ 2 j — L n jG n+ l t n-l — B n jG ntn -i, 71 > 2. 

It has some interest to remark here that, as described in |9j, since 

(72) rank(L nj ) = n + 1 = rank(C n+ i i:? ) , j = 1, 2, tj > 0, 
the columns of the joint matrices 

L n = (-^n.,i ) L n 2 ) and (7 n = (C/ n>1 , C n 2 ) 
of size (2n + 2) x (n + 2) and (2n + 2) x n respectively, are linearly independent, i.e. 

(73) rank(L„) = n + 2, rank(C n ) = n. 

Therefore, the matrix L n has full rank so that there exists an unique matrix D^ n of size 
(n + 2) x (2n + 2), called the generalized inverse of L n : 

(74) Di = (D nA \D n , 2 ) = (L T n L n y 1 Ll 
such that 

D n L n = I n +2- 
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Moreover, using the left inverse of the joint matrix L n 

(\ \ 

1/2 o 1/2 o 



O 



V 



1/2 




1/2 



we can write a recursive formula for the monic orthogonal polynomials 



(75) 



P 



n+l 



Dl 



x 

y. 



F n -DlC n F n ^, n>0, 



with the initial conditions P_i = 0, Po = 1, where £g> denotes the Kronecker product and 
(76) B n = (-Bji , B^ 2 ) , C n = (Cn >x , C^ 2 ) ) 

are matrices of size (2n + 2) x (n+l) and (2n + 2) x n, respectively. This recurrence (1751) 
gives another representation of [9, (3.2.10)], already presented in the bivariate discrete 
case in [281. 



Corollary 5.3 (Structure relations). For monic polynomials, n > 3 and j = 1,2, with 
the notation (f^j for d>j(x. y) given in (f^Ijj. i/ie coefficients of the structure relations 
are given in terms of the coefficients of the second order partial differential equation ( TJ1| ) 

(77) = E nj - (a^Ln-i^L^i + (5jL n ^i^L n) i + ^jL^i^Ln^) , 

(78) 5 nj - = E n j (5jL n _i i + £ J L n _i j2 ) — W n jG n+ i tn 

+ (j n) „_iE n _ij (ajL n _ 2)1 L n _ lil + /3jL„_2,2-^n-i,i + 7^-^0-2,2^-1,2) > 

(79) T nj - = t>JjE n j + G n) „_iE n _ij (5jL„_2,i + SjL n ^2,2) — W nj jG n +i tn -i — Sn,jG n ^ n —\ 

+ G n!n _2E n _2j (ajL n _3 j iL n „2,l + 0jL n -Z^L n -2,l + 7i-^n-3,2-^n-2,2) • 

Corollary 5.4 (Derivative representations). For monic polynomials, n>2 and j = 1,2, 
the coefficients of the derivative representations (E3|) are given in terms of the coefficients 
of the second order partial differential equation fOfy by 

(80) = (LnjEn+ij) -1 , 

(81) = (f?n,n-l ~~ Ki,j-^nj'^n+l,nEnj)Ki-lJj 

(82) Z n j = (G nin _2 — Ki,i-^n,i^n+l,n-lEn-lJ — ^y'^n-ljG'fi.n-l^n-lj) Ki-2J- 

The results obtained in Section H] and [5] can be applied to any polynomial solution of 
an admissible partial differential equation of hypergeometric type. In the next section 
we study in detail the monic Appell polynomials [I]. The non-monic case is also briefly 
discussed. 
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6. Example: Monic Appell polynomials 

In this section we give explicit expressions of all the aforementioned matrices appearing 
in the three-term recurrence relations and structure relations, by applying the results ob- 
tained in the previous sections to the case of monic Appell polynomials. Other examples 
of orthogonal polynomials of two variables can be treated in a similar way. 

In 1882, Appell jl] introduced a three parameter family of polynomials of degree n + m 
in terms of generalized Kampe de Feriet hypergeometric series [31] which in the monic 
case and reduced to two parameters can be written as [TT| Equation (12), p. 271] 

(83) A^(x,y) = (-1)^ J?[ nW h, 

(a + (3 + n + m) n+m 

Hi-vx ( a + (3 + n + m : -n; -m 

X *0:1;1 o X,y 

y -:a]P 
where a > and > and (a) n denotes the Pochhammer symbol. 

6.1. Second order partial differential equation. The admissible partial differential 
equation of hypergeometric type satisfied by A^£\x,y) is [TTJ Equation (15), p. 272] 

(84) ,(l - - 2 xy d -^ + ,(l - yf-j^- + {a -( a + P + 1),)^- 

+ [fi - (a + $ + l)y)—^- + ( n + m)(a + f3 + n + m)A%$ = 0. 

oy 

From (J23J) we obtain that the partial derivatives 

^H-,y) = ^A^)(x,,), 

satisfy the admissible partial differential equation of hypergeometric type 

d 2 d 2 d 2 

x(l - x)—z {r ' s \x,y) - 2xy-^^z {r ' s \x,y) + y(l - y)—z ir ' s) (x,y) + 

d 

+ (a + r - (a + (3 + 1 + 2(r + s))x)— z {r ' s \x, y) 

ox 

+ ((3 + s-(a + (3 + l + 2(r + s ))y)-^z^ s \x, y) 

oy 

+ (n + m — r — s)(a + f3 + n + m + r + s)z^ r ' s ^ (x, y) = 0. 

It is easy to check that, from the above differential equation, the matrices G nn ^i and 
G ntn -2 in the expansion 

(85) A n = A^(x, y) = (A^f(x, y), . . . , A^g(x, y), . . . , 7$f(x, y)) T 

= x™ + G n , n ^ix. n 1 + • ■ ■ + G n fl x° , 

are explicitly given from the general expressions fl65l) and fl66l) . 



18 AREA, GODOY, RONVEAUX, AND ZARZO 

6.2. Orthogonality. From f|T9|) we obtain the following weight function 

(86) e(ar,y) = x a -y- 1 . 

From fj!4H the positivity of g(x,y) yields the following triangular domain 

(87) K = {(x,y) : x > 0, y > 0, x + y < 1}. 
Therefore, the orthogonality property reads as 

g(x, y)A%£\x, y)k { «f (x, y) dx dy = A^8 N , K , 

n 

where N = n + m, K = k + s. 

6.3. Three term recurrence relations. For n > 0, monic Appell polynomials satisfy 
the three term recurrence relations 

xA n = L n) iA n+ i + B nj xA n + C nj iA n _i , 

yA n = L nj 2A„ + i + B nt 2A n + C n ^A n ^i , 

with the initial conditions Aq — 1 and A_i = 0, where A n is defined in ( 1851) and L n j are 
defined in (jUJ). Using fl68|) and fl69|) the recursion coefficients B n j are given by 



(89) 



where 



o,o 



bn,n—l b n ^ n J 



(n — i)(a + n — 1 — i) (n+ 1 — i)(a + n — i) 



and 



(90) 



with 



2n - 1 + a + (3 

2(i + l)(P + i) 

(2n - 1 + a + (3)(2n + 1 + a + fi) 



( b , 6o,i 
61,1 61,2 



2n + 1 + a + (3 

< i < n - 1 



< i < n, 



B 



n,2 



b n -l,n-l &n-l, 

bn..n 



bu. 



J i,i+1 



v O 

z(2n-z + a) (z + l)(a + 2n + 1 - z) 



2n-l + a + [3 2n + 1 + a + /3 

2(n — i)(a + n — 1 — i) 



0<i<n, 



(2n - 1 + a + /3)(2n + 1 + a + (3) ' 



< i < n - 1. 
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Moreover, using ( 1701) and fITTT) we have 
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(91) 



( c o,o 

c l,0 c l,l 

c 2,0 c 2,l C2,2 



o \ 



Cn— l,n— 3 ^n— l,n— 2 Cn— l,n— 1 
Cn,n—2 Cn,n—1 / 



where 



_ (n — i)(a + n — 1 — i)(n + i + f3)(n — 1 + i + a + f3) n ^-^ 



and 



(2n + a + /3)(2n - 1 + a + /3) 2 (2n - 2 + a + 0) 
(z + 1)(J3 + z)(2(n - z - l)(n + z + /?) + a{2n + a + - 2)) 

(2ra + a + 0){2n - 1 + a + /3) 2 (2n - 2 + a + /3) 

(z + 2)(z + l)(^ + z')(/3 + z + l) 



< z < n- 1, 



(2n + a + (5)(2n - 1 + a + /3) 2 (2n - 2 + a + 0) 



< z < n - 2, 



/ Co,0 c 0,l 
Cl,0 Ci,i 



(92) 



n,2 



Cl,2 



2 C n _2,n-1 
Cn— l,n— 2 Cn— l,n— 1 
Cn,n— 1 



with 



(n - z)(a + n - 1 - z)(/9(2n - 2 + /?) + a{2i + (3) + 2%{2n - 1 - z)) 
(2n + a + /3)(2ra - 1 + a + /3) 2 (2n - 2 + a + /3) ' 
(z + l)(a + 2n - 1 - z')(/3 + z')(a + /3 + 2n - 2 - z) 



(2n + a + (3){2n - 1 + a + f3) 2 (2n - 2 + a + f3) ' 
for < z < n — 1 and 

(n — z)(n — 1 — i)(a + n — 1 — z')(a + n — 2 — i) 
= (2n + a + /3)(2n - 1 + a + /3) 2 (2n - 2 + a + (3) 



for < z < n - 2. 



6.4. First structure relations. The partial differential equation (184)) for monic Appell 
polynomials corresponds to cases (vi), (ix) and (x) in Section |3~T1 Therefore, we obtain 
that 



(93) 



^ s \x,y) = [x(l -x- y)Y [y(l -x- y)] s . 
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As a consequence, the structure relations fl60|) satisfied by monic Appell polynomials 
defined in ( 1S5|) are given by 



(94) 
(95) 



x(l — x — y) 



d_ 

dx 
d 



W n i A n+ i + S^iA,,, 



y n—l: 



y(l-x- y)—K = W n}2 A n+1 + S n>2 A n + T n . 2 A n _i, 
dy 



for n > 1, where using f !77|) we get 



U>ii W h2 



(96) 



With Wj 



n + i = Wi < i < n, and 



/ wo,o ^o,i 
Wii w 1>2 



(97) 



with w,- 



% = Wa+i, < % < n. 



Moreover, from (1751) we obtain 



w n _i, n _i w n _i, n 

w n,n w n,n+l / 



( s 0,0 s 0,l 
Sin Sii 



(98) 



s l,2 



"Sn— l,n— 2 Sn—l,n—l ^n—l,n 

/ 



with 



H,i+1 



(n - z')(-n + (2n - l)i - 4z 2 + (n - 2 - 3i)/3 + a(n - 1 + i + a + (3)) 
(2n - 1 + a + (5){2n + 1 + a + (3) ' 
(n - i)(n + a)(2i + 1 + a + 0) 



(2n - 1 + a + 0)(2n + 1 + a + /3) 



for < z < n — 1, and 



2(i + l)(n - 1 - + i) 
(2n- l + a + /3)(2n+ l + a + /3) : 
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for < i < n — 2. Also, using f JTHj) we have 

/ Q \ 

Sl,0 s l,l $1,2 
S2,l S2,2 
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(99) 



>n,2 



V 



Sri— l,ra— 2 "Sn— l,n— 1 "Sn— 1,1 
Sn,n—1 Sn,n ) 



with 



i(P - (3 2 - i + (3 i + 4 i 2 - a (-2 + + 3 i - 2 n) - 2 (-1 + + 3 i) n + 2 n 2 



for 1 < z < n, and 



s i+l,i 
s i,j+l 



(2n + 1 + a + 0)(2n -1 + a + fi) 

(1+j) + + Q + 2_z + 2n) 

(-1 + a + /3 + 2n) (1 + a + /3 + 2ra) 
2 i (— i + n) (—1 + a — i + n) 



[-l + a + (3 + 2n) (1 + a + /3 + 2 n) ' 
for < i < n — 1. Furthermore, from (1791) it holds 



^0,0 


^0,1 


^1,0 




^2,0 


^2,1 



(100) 



^1,2 
^2,2 



L 2,3 



where 



^n-2,n-4 

£71— 1, 71— 3 "n— 2 ^n— 1,tj— 1 

v o 00 

(n — z") (n — 1 + a — z) 



*i+l,i 



(-2 + a + /3 + 2 n) (-1 + a + (3 + 2 nf (a + (3 + 2 n) 
x (/3 2 (1 + z) + z 2 (1 + 3 z) + a (3 (1 + n) + a 2 (-z + n) 
+/3 (z (3 + 4z) +n (-2z + n)) + n (- ((-2 + z) z) + n (-1 - z' + n)) 
+a (i (2 + z) + n (-1 - 2 i + 2 n))) , < i < n - 1, 

(n — z)(n — z — l)(a + n — 2 — z)(a + n — z — l)(a + /3 + n + z) 



(2n + a + /3)(2n - 1 + a + /3) 2 (2n - 2 + a + /3) 
(/9+i)(n-i-l)(i+l) 



0<i<n-2, 



U 



(2n + a + (3)(2n - 1 + a + /3) 2 (2n - 2 + a + /3) 
x (a(a + /3 + n + z - 1) + (3{n - 2z - 3) + (-2 + (2n - 5)z - 3z' 2 )), < z < n - 2, 
(0 + + i + l)(n - i - 2)(i + l)(i + 2) 



i+2,t 



(2n + a + (3)(2n - 1 + a + (3) 2 (2n - 2 + a + (3) 



< z < n - 3, 
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(101) 



where 





/o 


O ) 










to n 


to 1 to 9 £9 ^ 
z,i z,z z,o 


T n ,2 — 










^ri— 2,n— 4 ^n— z,n— 3 b n—2,n—2 z,n— 1 






^71 — 1,72 — 3 71 — 1,72 — 2 ^71 — — 1 




V 


(f^j tn,n—2 tn,n—l J 


i(- 


-i + n 


) (— 1 + a — i + n) 


(-2 + a + + 2n) (- 


1 + a + + 2nf (a + f3 + 2n) 


(a0 + 2 -i (l + 3i- 


-An) -0 (2 + i-2n) + a (-l + 2z- 


C. i < n — 1, 






2 (—1 — z + n) 


(— z + rz) (—2 + cc — z + n ) (—1 + a — z + n 



-2 + a + /3 + 2n) (-1 + a + /3 + 2n)^ (a + /3 + 2n) 
(1 + i) (/9 + i) 



1 < i < n - 2, 



H+2A 



(-2 + a + + 2 n) (-1 + a + + 2 n) (a + + 2 n) 
x (-3 (1 + if + (1 + n) (a + n) (a + + 2 n) + (1 + if (1 + 4 a + + 8 n) 
+ (a (3 + a) - (-2 + 0) + 4n + 6 a n + 6 n 2 )) , < i < n - 1, 
(1 + i) (2 + i) (0 + i) (1 + + i) (-2 + a + - i + 2 n) 



-2 + a + + 2 n) (-1 + a + + 2 nf (a + + 2 n) 



0<i<n-2. 



6.5. Derivative representations or second structure relations. The monic Appell 
polynomials defined in ( 18"5"|) satisfy the derivative representations 



(102) 

where Qn 



l n = V n>1 QW + y nil Q^! + Z n>1 Q^_ 2 , n>2, j = 1,2, 



L„,i4:A„+i and 



defined in ( 180|) are given by 



f v o,o 



(103) V n , 



and 



(104) K 



L n ,2^^n+i- In this example the matrices V Ui 
\ 



with v 



n + 1 — z 



< z < n. 



n,2 



via 



v n-l,n-l 

V r , 



with £>,- 



z+l' 



< i < n. 



BIVARIATE PARTIAL DIFFERENTIAL EQUATIONS AND OP SOLUTIONS 

Moreover, from flHTj) we have 
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(105) 



where 



A/o,o 
2/i,o 



O 



Vn-\,n-2 Vn-\,n-\ 
Vn,n-1 ) 



Vi,i 
Vi+i,i 



2i + 1 - a + (5 



(2n + 1 + a + P)(2n - 1 + a + 0) ' 

2(i + l)(/3 + i) 

(n - i)(2n + 1 + a + 0)(2n - 1 + a + 0) ' 



< i < n- 1, 



and 



(106) 



where 



/j/o 



n,2 







2/1,0 2/i,i 



V O 



\ 



2/n-2,n-3 Vn-2,n-2 

Vn-l,n-2 2/n-l,n-l 
2/n,n-l / 



3/i+M 



2(n — i)(n — 1 — i + a) 



(1 + i)(2n + 1 + a + /3)(2n - 1 + a + /3) 
2n-l-2i + a-(3 



(2n + 1 + a + P)(2n - 1 + a + 0) 



< % < n- 1. 



Also, from (EH) we have 



(107) 



Ai,l 



^0,0 

^1,0 #1,1 
^2,0 Z 2A 



z 2,2 



z n~2,n-4 z n-2,n-3 ^n-2,n-2 
z n— l,n— 3 ^n— l,n— 2 
^n,n-2 y 
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where for < i < n — 2, 

(n — i)(n — 1 — i + a){n + % + (3) 
Zi,i ~ ~ (2n + a + (3) (2n - 1 + a + (3) 2 (2n - 2 + a + 0) ' 
_ (i + l)(-2(i + 1) + a - P)(P + i) 

Zi+1,i ~ (2n + a + P)(2n - 1 + a + P) 2 (2n - 2 + a + (3) ' 

(z + l)(z + 2)(/3 + z)(/3 + z + l) 

^ +2,i ~ (n- 1 -i){2n + a + P){2n- 1 + a + /3) 2 (2n - 2 + a + /3) ' 



and 



' z ,o 

z 2,0 



(108) 



^n,2 



\ 



Zl,l 
%2A 



z 2.2 



Z n~2,n-A Zn-2,n-3 z n-2,n-2 
z n— l.n— 3 *n— l.Ti— 2 







^n,n— 2 / 



where for < % < n — 2, 



(—1 — i + n) (— i + n) (—2 + a — z + n) (—1 + a — i + n) 
(1 + i) (-2 + ct + /3 + 2n) (-l + a + /3 + 2n) 2 (ct + /3 + 2n 
'(-1 - i + n) (-2 + a - i + n) (a- (5 + 2 (-1-i + n)) 
(-2 + a + (3 + 2 n) (-1 + a + (3 + 2 nf (a + P + 2 n) 

(2 + i) (1 + + g (-2 + a - i + 2 ra) 

'-2 + a + £ + 2 n) (-1 + a + P + 2nf (a + P + 2n)' 



6.6. Non— monic orthogonal solutions of ([84]). For any non-monic orthogonal poly- 
nomial solution of the partial differential equation (IMj) it is possible to obtain the main 
differential and algebraic properties by using the results given in Section HJ In this Sec- 
tion we give some relations for two concrete non-monic solutions of ( 18"4"|) . orthogonal with 
respect to fl8"B"|) in the domain (J8"7|) . 

On one hand, also in 1882, Appell considered a family of non-monic polynomials 
solution of the partial differential equation (j8"4"|) . This orthogonal family can obtained 
from the Rodrigues formula (jUJ) (see [TH Equation (11), p. 271]) using the weight (|93|) 



(109) F^\x,y) 



x 1 a y 1 ~ 



gal 



[x 



n+a-1 y m + 



-l(l_ x _ y )"+m] . 



(<*)„ (P) m dx n dy* 

These polynomials can also be obtained from the classical Appell's orthogonal polyno- 
mials defined in [32], Equation (6), p. 63] by taking 7 = a + p. 

Clearly, both, monic Appell polynomials defined in fl83|) and the non-monic family 
(11091) . form a biorthogonal system in the domain flH7|) with respect to the weight function 
(USD pU Equation (17), p. 272]) 



I x^y-'F^ (x, y) Aff> (x, y) dx dy = 5 nk 5 ml A 



a,P>0. 
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Then, if we denote 
(110) F n = Wf»»{x,y) = (F { :f(x,y), . . .,F^{x,y), . . . , F%f{z, y)f 

we have the following formula linking both solutions in column polynomial vector form 
(HI) F n = G£ n A n , 

where A n is defined in (1531) and the entries of the matrix G^ n = (gfj(n)) of size (n + 
1) x {n + 1) have the following explicit form 

ft >> = U — — • - hl - 

Once the matrix is known, applying the formulae given in Section H] it is possible 
to obtain the coefficients of the thee-term recurrence relations, structure relations and 
derivative representations for this non-monic family (I109p . 

On the other hand, let us consider the non-monic Koornwinder family defined in [91 
Section 2.4.2, p. 86] with a -»■ a - 1/2, -»■ (3 - 1/2 and 7 = 1/2: 

(112) K$P(x,y) = P^+^)(2x-l)(l-xrP^- 1 \^ - 1), > 0. 

As already mentioned, this family is another polynomial solution of the partial differential 
equation ( 1541) orthogonal with respect to ( 156]) in the domain ( 157|) . where Pn (x) are 
the Jacobi polynomials [H]. Therefore, monic and non-monic Appell polynomials and 
this family form also biorthogonal systems. If we denote 

K n = K^\x, y) = (K^f(x, y),..., K^g(x, y), . . . , J$f(x, y)) T 
then we have 



(113) K„, = G« n A n , 

where the entries of the matrix G^ n = (g^j(n)) of size (n + 1) X (n + 1) have the following 
explicit form 

(0, i<j, 
ftjM = < (a + P + n + i^-iiP + fii . > . 
( (n - i)\ j\ (i - j)\ 

for < i, j < n. 

Relations fillip and ( II 13[) between monic (153]) and non-monic (I109P Appell polyno- 
mials, and between monic Appell ( 18 3 j) and Koornwinder (I112p polynomials solve the 
following connection problems 

n 

F^Jlix, y) = V g?An)AL a l®(x, y), < I < n, 



j=0 
n 



K^ e (x,y) = £<£(n)X£f>,y), 0<£<n, 

3=0 

between the corresponding scalar orthogonal polynomial families. 



26 



AREA, GODOY, RONVEAUX, AND ZARZO 



Finally, we would like to mention here that our goal is not to exploit all possible 
situations covered by our approach, but to emphasize its systematic character, which 
allow one to implement it in any computer algebra system, here the Mathematica [33] 
symbolic language has been used. 
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